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The Talmi-Moshinsky transformation is applied to the strong decay and production of mesons. 
Amplitudes for conventional mesons are obtained in a simple way and are generalised to encompass 
the production and decay of hybrid states. The selection rule disfavouring hybrid decay into pairs 
of identical S-wave mesons has a natural explanation and a new class of selection rules is uncovered 
with implications for production of the exotic 1 _+ hybrid in e + e~ experiments, charmonia decays 
and the lattice. 
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1 I. INTRODUCTION 

O ' 

£S) , The flux tube model sets a scale for hybrid decays relative to those of conventional mesons 0, H, H[ , a process which 
• the lattice is now able to calculate The agreement between the two is good in a simple version of the model 
[ with equal wavefunction widths and a simplified string overlap term [f|. A decay formalism is developed here which 
works within the same level of approximation and treats the production and decay of states with arbitrary gluonic 
excitation. A technique borrowed from nuclear physics - the Talmi-Moshinksy transformation - provides a simple 
framework for the calculation and uncovers new selection rules sensitive only to the gluonic angular momentum of 
the initial state. Using the effective quark-flux tube wavefunctions developed elsewhere Q decays involving hybrid 
, states are a straightforward extension of those involving only conventional mesons. 

Within this formalism the amplitude for the decay of any meson into any other pair of mesons is obtained, be they 
£NJ , orbitally, radially or gluonically excited (hybrid) states. Previous results for conventional and hybrid meson decays 
C*~) ■ are derived in a simple way and new results obtained for hybrids in the final state. The amplitudes are left in rather 
general form: simpler polynomial forms will be collated elsewhere and discussed with reference to the phenomenology 
of hybrids in experiment and on the lattice Q ■ The selection rule against the decay of a hybrid meson into pairs of 
identical S-wave states can be generalised using the symmetry of the Talmi-Moshinsky transformation, giving several 
, new results. A new selection rule forbids the decay of any meson with the gluon field in its ground state to the 
1 ^ hybrid and a IS meson in the limit of equal width outgoing meson wavefunctions. Within model assumptions 
developed elsewhere this rule may help discriminate interpretations of the A" (3940) and F(4260) states, and modes 
such as ip — > wip and x ~~ > ^l 71 " are forbidden. 
l The Talmi-Moshinksy transformation and its relevance to quark-pair creation processes is discussed in section |TT] 
and applied to the decay of a conventional meson to any other pair of conventional mesons in section [TTT1 Section HVl 
introduces effective wavefunctions for mesons with flux tubes which and these are applied to the decay and production 
of hybrids in sections fVl and PVTl In section TVIII some phenomenology of hybrid production in e + e~ experiments is 
discussed. 
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II. QUARK-PAIR CREATION AND THE TALMI-MOSHINSKY TRANSFORMATION 



Central forces in nuclear physics require matrix elements of the form 

((mh <g> n 2 k)j 12 m 12 |^(|ri - r 2 |)|(ni/i ® n' 2 l' 2 ) v ^ m Q (1) 
where the coupled wavefunctions in the nucleon coordinates ri , r 2 are 

(ri,r 2 |(ni/i <g> n, 2 Z 2 )i 12 m 12 ) = ^ (ri|niZimi)(r 2 |n 2 Z 2 m 2 )(Zimi , / 2 m 2 |Zi 2 mi 2 ). (2) 

mim 2 

In the basis of equal width harmonic oscillator wavefunctions there is a convenient procedure for the evaluation of 
such matrix elements. The wavefuctions in (ri,r 2 | are expanded in a finite sum of wavefunctions in vectors (R, r| 
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FIG. 1: Strong decay topology. 



proportional to the relative and centre of mass coordinates of the nucleons, 
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V2 ' 
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r i + r 2 

V2 ■ 



(3) 



where to match onto later notation the conventions of upper and lower case to denote the relative and centre of mass 
coordinates are the opposite of those normally used. Denoting by upper case letters the quantum numbers in R and 
lower case letters the quantum numbers in r, the expansion gives 

{{nil 1 ®n 2 l2)i 12 m 12 \V{\Yi -V2\)\{n[l[®n' 2 l' 2 ) Vi2m , i2 ) = 

Y, ^ nih ® n 2h)h 2 \(NL <g> nl) h2 )((NL ® nl) ll27ni2 \V(V2R)\(N'L' <g> n'l'^^J 



NLnl N'L'n'V 

((N'L'® n n') l , 2 \(n' 1 l[®n 2 l 2 ) Vi 
The centre of mass co-ordinates integrate out entirely leaving only a simple one dimensional radial integral, 

((NL <g> nl) h2mi2 \V{V2R)\(N'L' ® n'l')i' l2 m' 12 
The coefficients in the coordinate transform, 

,ri - r 2 ri + r 2 



(4) 
(5) 

\(NL®nl) h2mi2 ) 
(6) 

are the Talmi-Moshinsky brackets [g, Q whose values are well known [l(| • The essential idea is that the harmonic 
oscillator Hamiltonian can be written equivalently in either pair of Jacobi coordinates (ri,r 2 ) or (r, R): each con- 
struction shares common eigenvalues in li 2 2 , and H. The latter demands that the number of energy quanta is 
conserved, 



((nth ® n 2 l 2 ) h2 \(N L ® nl) ho ) = 1 V" / d 3 rid 3 r 2 ((niZi <g> n 2 / 2 ) /l9mi , |n, r 2 ) (- 

ll\ 2 + l '-^ J 



V2 ' y/2 



2ni + h + 2n 2 + h = 2N + L + 2n + I, 



(7) 



so the summations in ^ are always finite. 

In quark pair-creation processes the dominant topology is that of FIG.[T] Hairpin diagrams are generally suppressed 
or forbidden on account of colour considerations, or in the flux tube model because the created pair sit at the ends of 
two broken pieces of string constituting the final state hadrons. The corresponding geometry is show in FIG. [2] where 
r,ri,r 2 are the qq coordinates of the initial and final states respectively, and 



R 



1'2 



r = ri + r 2 . 



(8) 



The spatial amplitude for such processes has integrals over the wavefunctions of the initial and final states with 
coordinates r,ri,r 2 and the relative wavefunction of the outgoing mesons in the relative coordinate r/2. More 
generally, allowing for a transverse degree of freedom R (such as in a hybrid meson) the integrals are of the form 



J d 3 r!d 3 r 2 d 3 rd 3 R(5 3 (r - n - r 2 )<5 3 (R - 



ri — r 2 , 



(9) 




FIG. 2: Strong decay geometry. 
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Only two of the integrations need to be performed manually on account of the delta function: the only complication 
that arises is that at least one of the wavefunctions has to be expressed in the coordinates of another. This can be 
done using bipolar expansions but the presence of operators connecting the initial and final states, particularly in 
cases involving radially excited wavefunctions, makes the algebra unwieldy and masks underlying symmetries. A more 
convenient approach is to note that in its nature the problem is essentially the same as the earlier nuclear physics 
one: the wavefunctions must be translated from one set of Jacobi coordinates to another. Apart from a scale factor of 
s/2 the pair of coordinates |(5J) is that obtained from the orthogonal coordinate transformation on the vectors (ri, r 2 ) 
cf. ([5]) above. Using the operator formalism of the harmonic oscillator the matrix elements for such processes can be 
evaluated without having to perform any integration manually apart from simple radial integrals at the final stage. 

The approach is to first operate on the wavefunctions of the initial or final state, leaving wavefunctions in the 
coordinates r 1; r 2 and r (and in more general cases R); those in the former coordinates are translated via the Talmi- 
Moshinsky transformation into a new set of wavefunctions of the latter coordinates, 

((mh ® n 2 l 2 ) ll27ni2 \r u r 2 ) = ^ {{mh <8> n 2 l 2 )i 12 \{NL <g> nl)i 12 ){(NL ® nl)i 12mi2 \V2R, r/V2), (10) 

NLnl 

leaving only simple radial integrals. There are a finite number of ways of allocating the degrees of freedom to the 
transformed state, owing to the energy conservation rule ([7]) and its corollary, 

(_)£+' = (_)'H-'a. (11) 

This approach is well suited to the problem of conventional and hybrid meson decay: the translated wavefunction has 
a unit of orbital angular momentum L in the transverse degree of freedom R, which must correspond to the initial 
transverse orbital angular momentum - zero for a conventional meson or one for the first excited hybrid state. Along 
with the symmetries of the Talmi coefficients, 

{{ruh ® n 2 l 2 ) h2 \{NL g> nl) h2 ) = (-) I+I ' J {{n 2 l 2 g> ni h)i 12 \{NL ® nl) ll2 ), (12) 

this uncovers a new class of selection rules that probe the transverse degrees of freedom of the initial state, of which 
the simplest example is the suppression of hybrid decays to pairs of identical S-wave mesons. 

Most decay models are nonrelativistic in character and assume L, S factorise: the pictures differ in the assumed 
quantum numbers of the created pair. An early model assigned the emergent pair the quantum numbers of the 
vaccum, corresponding to qq pair creation in a 3 Po state [ll|. With the arrival of QCD it was proposed that instead 
the qq should have the 3 Si quantum numbers of the vector gluon but this was found to be less successful in describing 
the data [12]. In the strong coupling expansion of lattice QCD the natural degrees of freedom arc not quarks and 
gluons but quarks and flux tubes and it is not a pcrturbative gluon that triggers decay but the breaking of a gluonic 
flux tube [l[. The Hamitonian allows for the quarks at the end of the broken pieces of string to carry either 3 P or 
3 Si quantum numbers with corresponding operators: 

3 Po : Xi ' V = xi ■ (V ri + V r2 + V r/2 ) (13) 
3 S i: xi -f (14) 

where %i denotes the spin 1 wavefunction of the emergent pair and the vectors r, ri and r 2 are those of FIG.O In the 
flux tube model of Isgur and Paton it is argued that the zero point oscillations of the string wash out the 3 Si operator 
and that the 3 Po operator survives, converging on the phenomenologically successful 3 Po model and embedding it in 
a more appealing physical picture [13j . An additional feature of the flux tube model is that the amplitude for decay is 
driven not only by the overlap of the quark wavefunctions but also by those of the flux tubes of each hadron. In the 
next section the amplitude for the 3 Po decay of any conventional meson into any other pair of conventional mesons 
will be formulated, neglecting the string degrees of freedom. For calculations involving hybrids the string degrees of 
freedom must be included: this is presented in the subsequent sections. 

III. CONVENTIONAL MESON DECAY 

The relevant degrees of freedom are the colour, flavour, spin, orbital, radial and total angular momentum quantum 
numbers of each meson. The colour factor is the same for all processes of the type in FIG. [1] and contributes only an 
overall factor which is absorbed into the normalisation. Flavour overlaps can be handled by standard techniques and 
will not be discussed here. For the initial meson with its qq axis at r, the wavefunction assumes the form 



(r|(s ® nl)jm) 



(15) 
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where s corresponds to a pair of quarks whose labels have been suppressed coupled to spin s and the space part is 
given by 

(r\nlrrn) = (r|n0^ m '(f)- (16) 

The radial wavefunctions are taken as harmonic oscillators governed by a width parameter /3, which in the present 
work is taken to be universal; typically j3 ~ 0.4GeV for light quarks. Corrections for wavefunctions of unequal 
width can be implemented at a later stage by mapping more realistic wavefunctions onto linear combinations of these 
harmonic oscillator basis states. Taking a universal wavefunction width suggests a system of natural units in which 
P = 1, so that the normalised radial wavefunction is 



With the quantum numbers as per FIG. [1] and in the rest frame of the initial state, the final state consists of a pair 
of mesons ((si ® mZi^-Jri) and ((s 2 ® n 2 l 2 )j 2 l r 2) with relative coordinate r/2 moving apart with (dimensionless) 
momenta of magnitude p measured in units of the wavefunction width (3. The final states are coupled to ji2 and are 
in a relative partial wave L, 

(LM L \r/2) = (r/2\L)Y^(r), (18) 
where in the absence of final state interactions the relative radial wavefunction has the form 

(L\r/2) = ^i L j L (pr/2). (19) 

The final state must be coupled to angular momentum \jm) on account of its being connected to the initial state by 
a scalar operator, so has the form 

((((si ®nih) h ® {S2®n 2 l2)j 2 )j 12 ® L) jm \n,T2, r/2) (20) 

The total width for the decay of a conventional meson into another pair of conventional mesons is the sum over 
couplings ji2 and relative partial waves L of the squared amplitudes: 

T(snlj -> swihji + s 2 n 2 l 2 j 2 ) = ^ll^jj^- ^{{{{s\ ® n ^)h ® ( s 2 ® n2h) j 3 )h 3 ® L )j\\xi • V||(s ® nl)j) 2 (21) 

jviL 

where here 70 is a phcnomenological pair creation constant. Details of choice of phase space will not be discussed in 
this paper, which is concerned only with the evaluation of the matrix element itself. In the matrix element above, 
and in subsequent manipulations, the state vectors |ri,r2,r/2) and (r| have been dropped and an overall integration 

J d 3 nd 3 r 2 d 3 rd 3 R(5 3 (r - n - r 2 )(5 3 (R - Tl ~ * 2 ) . . . (22) 

is implied. The first step is to separate the spin and space degrees of freedom; the relevant recoupling is essentially 



the same as that of 14| but it is useful to manipulate the radial quantum numbers along with the orbital parts, 



(((0i ® n x h) Sl ® (s 2 ® n 2 l 2 ) n ) n2 ® L^Hxi ■ V||(s ® nl) 3 ) = ^ (-) s+L+Sl2+ ' 12+i/ n ;/sl2;i2JU2J12 

512^12^/ 



si h ji 
s 2 h j 

Sl2 h 2 jl2 



S2 h 32 U Sl L 2 h ? T f \{j i f )}((si®S2) s J\xi\\s)(((n 1 h®n2l2)i 12 ®L) lf \\nnl) (23) 



with 



Hob... = V / (2a+l)(26+l)... (24) 



The spin part is a 9-j coefficient, 



1 S S12 
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The space part can be evaluated as follows: 

(((mil ® nafejiu ®i)i / ||V||n/) = (((Mi ® ^2)112 ® i)i/l|V||((ni/i ® n^) Jia ® L')i)j^ ( 26 ) 

x £<(ni*i ® n 2 Z^J(iV0 ® n'Ji a ), ia ) (27) 
JVn' 

x<((JV0 ® n'iia)^ ® iOillnO- ( 28 ) 

In the first step (|26p . the operator V acts on the final state wavefunctions to produce a new set of states which, to 
match the initial state ket, must be coupled to I. The operator V is the sum of three terms V ri , V r2 , V r / 2 which act 
on the (nili\, (712/2! an d (L\ parts of the wavefunction respectively. The matrix element in (|26p is found by recoupling 
to sandwich the operators between the relevant parts of the wavefunction: 

(((nih®n 2 l 2 ) h2 ®L) lf \\n((n[l[®n'AK t £}w|V||2/) 

+ Sn£ L L '(-) h+l * +h2 { l l? I ^}(M 2 ||V||n^)). (29) 

The matrix element (L||V||L') can be found by operating on (L\ with the definition (fT8|) : 

(£||V||L-1) = ipVL, (30) 
(i||V||i+l) = -ipVL + 1. (31) 
The matrix elements ('ii7i||V||n' 1 Zj) and (712^2 llVHri,^) follow from those of the creation and annihilation operators, 

a=-^(p-ir), af = -^=(p + ir), (32) 

so that V = ip raises or lowers the number of quanta by one [IH, [l6| ■ The matrix elements of r are also recorded here 
for later use, 



(n/||V||nJ + l)= (nf||r||nf+l) = + l)(n + I + 3/2), (33) 

(nJ||V||n + lZ- 1) = {nl\\r\\n + 11 - 1) = -y/l(n+l), (34) 

and 



(nT||V||n/) = (n/||V||nT) (35) 
(nY||r||ni) = -{nl\\r\\n'l'}. (36) 



In the second step (f27|) . the harmonic oscillator states in |ri, r 2 ) are translated into states in |-\/2R, •*/ V2), as in EQN. 

{{n{l[®n' 2 l' 2 ) Vi2m ,Jr u r 2 )=Y,{^ (37) 

Nn' 

Here the orbital angular momentum in the harmonic oscillator state in a/2R, is forced to be zero on account of the 
initial state \nl) having zero orbital angular momentum in that direction; this leaves orbital angular momentum 1' 12 
in the direction r. What remains is a sum over radial quantum numbers n', N restricted by the energy conservation 
rule ([7]). This takes care of the integral 

J d 3 nd 3 r 2 <5 3 (r - n - r 2 )<5 3 (R - ^— -) (38) 
leaving only an integral over d 3 rd 3 R, represented by the final line |(25|) ; the integral over dR is trivially one, so that 

poo poo 

(((MD ® rc/Zi 2 V i2 ®L')i\\nl) = (l' 12 L'\\l) / dRR 2 (N0\V2R) drr 2 (n' l' 12 \r / V2) (L' \r / 2) (r\nl) (39) 
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with 

{l[ 2 L'\\I) = n, iaV <li 2 0,L'0|IO) (40) 

/ dRR 2 (N0\V2R) = 1— a N0K (l + 2K)\\ 

"'° K=0 



(41) 



where the index a = L' + / + 2A; + Z' 12 + 2k' + 3 and coefficients from the expansion of the Laguerre polynomials are 

R fc I> + / + 3/2) 
a„ ifc =^„/ (n _ fc)!r(fc + z + 3/2) - ( 43 ) 

The resulting expressions have a centrifugal barrier suppressing higher partial waves at small momentum and a form 
factor, 

p L x (polynomial in p 2 ) x e -15 " . (44) 

The polynomial form can be generated by repeated application of the confluent hypergeometric function recursion 
relation 

2 2 2 2 

M(a, b, - V -) =M(a-l,b, - ^M(o, 6 + 1, (45) 

and using A4(a, a, — fj) = e 15 5 this is always possible because the first two indices in M. will differ by an integer as 
([4"0")l demands 1' 12 + L' + I is even. The expressions so obtained recover the results of REFS. [l3| and [TtJ and extend 
upon them to consider decays not treated there, including modes that are not relevant the decays of physical mesons 
but are useful in modelling virtual qq decays 0] ■ 



IV. MESONS WITH FLUX TUBES 



In the flux tube model the state of the gluon field is manifest: a meson with angular momentum I is built from the 
quark (/') and flux tube (l s ) angular momentum 

l = l' + l s . (46) 

Because of the axial symmetry of the flux tube part of the wavefunction, the full wavefunction is not in general an 
eigenstate of l' 2 or l 2 and there are no associated quantum numbers l',l s . On the other hand, the projection of the 
angular momentum along the molecular axis is conserved quantity: for a state with qq axis at f the relevant operator 
is 1 • f = l s • f and the associated eigenvalue is A. The states with nonzero A are gluonic hybrids and the ±A states are 
degenerate left- and right-moving phonons, the interference of which gives parity doublets. The spatial wavefunctions 
can be described by kets of the form 

\nAl p mi ) (47) 

where here n is the quark radial quantum number, A = |A|, I is the total orbital angular momentum restricted to 
I > A, and P is a parity label relevant only to the A ^ states. In the parlance of molecular physics, states with 
A = 0, 1, 2 . . . are labeled E, II, A . . . in analogy with S,P,D. . . , the latter being used to label the quantum number I. 
The E states are the conventional quark model states; combining with spin gives the J PC quantum numbers: 

(48) 
(49) 
(50) 

The lowest lying orbital configurations for hybrid states are the UP^ states, with a single phonon in the lowest mode 
of excitation of the flux tube and total orbital angular momentum 1 = 1. Combining with spin this set includes the 
manifestly exotic 0" 1 , 1 ^ , 2^ states: 

| ?1 1 IIP+) 1++ |n 3 IIP£ 1)2 ) 0,1,2+- (51) 

In^PD I" |n 3 IIPo, 2 ) 0,1,2-+ (52) 



I^ESq) 


o-+ 


|n 3 ES!) 


1~ 




I^EPx) 


1+- 


l« 3 SP ,i,2> 


0,1,2++ 






2-+ 


l™ 3 £Di,2, 3 } 


1,2,3" 


he. 
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Those of most phenomenological interest are the quark radial ground states denoted |1IIP ± ) where here the prefix 
"1" stands for n + 1 as in the usual chemical notation. These are the family of hybrids which for light quarks are 
expected to lie at around 1.8-2.0GeV and are discussed in REFS. [3, 0|. As well as radial excitations |2ITP ± )&c., 
there exist in principle orbital excitations such as those with 1 = 2, 

|n 1 nD+) 2++ |n 3 nD + 23 ) 1,2,3+" (53) 

l^nrtf) 2— |« 3 nDr, 2 , 3 ) 1,2,3-+ (54) 

and then again a sequence of states with two phonons in the flux-tube |nAD+), | rt AF+)&c.[2l[ 

With the flux-tube degrees of freedom manifest the amplitude for a decay depends not only on the overlap of 
the quark wavefunctions but also on those of the flux tubes of the initial and final hadrons. Dowrick, Paton and 
Perantonis[l8j obtain the flux tube overlap for decays of the type 

E -► E + E, and (55) 
II -> E + E (56) 

for a harmonic flux tube connecting a fixed quark and antiquark. These factors provide a parameter-free scale for 
processes involving hybrids relative to those involving only conventional mesons, driven by calculable terms k and the 
ratio yfdj ' j3 where a is the string tension. For a wide range of other decay processes, including 

E n + E, (57) 
E -> n + n, and (58) 
A ^ E + E (59) 

the string overlap factors share the same generic formfi^: string breaking is suppressed exponentially according to 
the perpendicular distance p away from the qq axis, localising pair creation to a region of order and there are 

solid harmonic factors in p associated with the gluonic excitations II, A, &c. 

In the adiabatic approach the string overlap factors are used as a spatial weight for the full decay amplitude 0, [l3| . 
The subsequent calculations in this paper, which follow the prescription of REF. work in the approximation that 
the exponential localisation of string breaking is discarded. In this limit decay amplitudes involving only E states 
are those of the naive 3 Po model in which the amplitude for pair creation is equal everywhere in space; this is known 
to be a good approximation, essentially because the quark wavefunctions mimick this asymptotic form [13]. This is 
the level of approximation within which the amplitudes presented in Q were calculated, and the agreement with the 
lattice was found to be good. 

What remains of each string overlap is a solid harmonic associated with each gluonic excitation: in this sense 
the string overlaps are factorisable. In the approach of REF. [|| the appropriate factor is combined with the quark 
wavefunction of REF. [l[ to give an effective quark-flux tube wavefunction for each state; the overall decay amplitude 
is then driven by the factors k which are also calculated there. It is found that the combined quark-flux tube 
wavefunction can be expressed as a linear combination of tensors of rank I, 

(p,r|(A®n'0, ro ,> = 5^(p|AM)(r|n'l / m / >(AAi I iW|/mj>. (60) 

Jim' 

where the qq degree of freedom is a harmonic oscillator as before 

(r\n'l'm') = {r\n'l')Y^' (f) (61) 

and the flux-tube part of the wavefunction is 



(p|A ^ > = V2ATT pA ^ i(i5) - (62) 

The translation between the full wavefunction and the effective wavefunction, 

\nM p mi) ^^|(A®nY) ;mi )((A®nY) ; |nAZ p ), (63) 

n'l' 

is discussed in detail in REF. [||: here only the final results are quoted for the cases of most interest. For conventional 
mesons (A = E) the translation is trivial, 

\nlllmi) -> |(E®ni)i mi ) (64) 
= \nlmi) (65) 
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where in the second step, which follows from the definition (|60[) . there appears a harmonic oscillator kct defined in 
EQN. (Tl7o|) . Amplitudes involving only conventional mesons are then the amplitudes of the old 3 Po model with no 
string degrees of freedom, corresponding to pair creation with equal amplitude everywhere in space. 

For hybrid states (A = II, A, . . .) the effective wavefunctions are, unlike the full quark-flux tube wavefunctions, 
eigenstates of both V 2 and 1^ : 

I' 2 |(A®nY)* mi ) = l'(l' + l)\{A®n'l') lmi ), (66) 
l2l(A®n'0, m| ) = A(A + l)|(A®n'0, TOl ). (67) 

This is a consequence of having discarded the axially symmetric exponential localisation of string breaking. Thus 
£, IT, A . . . states have string angular momentum equal to the magnitude of its projection along the molecular axis, 
l s = A = 0, 1, 2, . . . The quark orbital angular momentum I' can take values |A — l\ < I' < A + / and the parity of a 
state is (— ) A+ ' +1 , so that the parity doublets correspond to /' even/odd. The quark radial wavefunctions for hybrids 
are subject to a modified centrifugal barrier, but decay calculations appear to be fairly insensitive to this [3|. In 
the approximation that the radial wavefunction is a IP harmonic oscillator state, the effective wavefunctions for the 
lowest lying hybrids are found to be 



inp-) -> |(II<8>lP)i) (68) 




|mp + > - ^/l^vlKn^isW + .-J+VH^V^ 1 ^ 01 ^^'-) (69) 

The ellipsis indicate contributions from higher radial excitations in the quantum number n' , owing to the mismatch 
between the orbital and radial characters of the quark part of the wavefunctions. For the orbitally excited states it is 
those with negative parity that pick up correction terms for the radial wavefunctions, which is assumed to be ID: 

|1IID-) «/r ( r./A|(n®lP) a ) + ...| + A /r ( ^ A /JL|(n®lF) a ) + ...| (70) 





5ir 

|1IID + ) |(II®1D) 2 ) (71) 

In a similar way it is possible to write down wavefunctions for A-type hybrids and so forth. 
Working in this basis, the decay amplitude that appears in equation (|21[) has the form 

((((si ® (Ai ® ni/i)ii)ii ® ( s 2 ® (A 2 ® n' 2 l' 2 ) h ) h ) jl2 ® L) 3 -||xi ■ V||(s® (A® n'l')i)j). (72) 

The spatial parts of the meson wavefunctions are still tensors of rank l,h,fo so the recoupling in equation ([23]) to 
separate the spin and space degrees of freedom is the same and the task is to evaluate matrix elements of the form 

((((Ai ® nili)ii ® (A a ® n' 2 l' 2 ) h ) h2 ® L) // ||V||(A ® n'l')i). (73) 

In formulating the full decay amplitude it is useful to associate p with one of the vectors appearing in FIG. O As 
discussed elsewhere @, the purely transverse oscillations of the string allows the replacement 

p -> R (74) 

for the flux tube vector of the initial state whose qq axis is at r. Within the harmonic approximation the same 
replacement can be made for flux tube excitations in the final state, and because the flux tube has only transverse 
oscillations this is equivalent to 

P - r 2 /2, (75) 
p n/2, (76) 

for the final state mesons whose qq axes are at ri , r 2 respectively. The quantum numbers and dynamical variables 
associated with each state are summarised in FIG. [3] The expansion (|60|) need never be done explicitly, and it is more 
useful to work in the coupled form (p, r|(A ® n'l')i m ). The solid harmonic (f6"2"| describes the elements of a tensor of 
rank A, so that for 

S, IT, A... (77) 

states the tensor A corresponds to 

1, p, {p®p} 2) ... (78) 
Thus apart from a scale factor the tensors associated with each gluonic excitation are: 
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(R,r|(A®nT) lm ) 




(r 2 ,ri|(Ai ® n' 1 l' 1 )i imi ) 



'T.m-ii 



(ri J r 2 |(A 2 ®n' 2 l 2 )i 
FIG. 3: State vectors and quantum numbers. 
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For the final states Ai, A 2 correspond respectively to tensors in r 2 , n which can be treated as operators on the \n' 2 l 2 ) 
and \nil'i) harmonic oscillator degrees of freedom that appear in (f?3")k this case is discussed in the next but one section. 
For an initial state hybrid A excitations are treated as operators on a harmonic oscillator degree of freedom in R 
that appears in the Talmi-Moshinsky transformation; this is discussed in the next section. For E states the matrix 
elements are just those of the unit operator, 

(ni||E||nV}=*nn'*H'nj. (79) 
For the II states the requisite matrix elements are those of the vector p, already given in (|33)) - (f36|) : 



H|II||nZ + l) = -v/(Z + l)(n + Z + 3/2), (80) 
(nZ||n||n+lZ-l) = -y/l(n+l), and (81) 
(nY||II||nZ) = -<nZ||II||nY). (82) 

For A > 1 hybrids the reduced matrix elements of the corresponding solid harmonic are required; these are known 
but such cases will not be discussed any further here. 



V. HYBRID MESON DECAY 



In this section the spatial part of the amplitude for the decay of a state of arbitrary gluonic angular momentum 
into a pair of conventional mesons 

A^S + S, (83) 

will be formulated, recovering the results of the previous section if the initial meson has A = S. The amplitude can 
be expressed in a form analagous to that of the conventional mesons case ([26]) - ([28)) : 

(((Mi <8> n 2 l 2 ) ll2 ) ® L) Z/ ||V||(A <S> n'l')i) = J2 ® "a^k, ® ||V||(«Zi g> n' 2 l' 2 \ 2 ® L')i>^ 4 ) 

x ((<l[®n' 2 l 2 ) l >J(NA®n"l"X 2 ) (85) 

Nn"l" 

x (((JVA ® n"l'\ 2 ®L)i\\(K® n'Oi). (86) 

Here the first step (|84|) is the same as for the conventional meson case (f26|) . In the second step the harmonic oscillator 
states in |ri, r 2 ) are translated into states in Iv^R, r /V^), as before, but now the translated states must have angular 
momentum A in the direction R to match the initial state |(A (g> n'l')i): 

((nili®^), iami jn,r a ) = Yl (K/i®^ 2 ^ 2 )r 1 J(^A®n'r) ri2 )((^A®n'r0r 12mi jV2R,r/V2). (87) 

Nn"l" 
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The matrix element in the final step can be integrated directly, but it is more convenient to treat A as an operator 
on the R degree of freedom and make use of the integrations already performed, 

(((NA®nin'%. 2 ®L) l \\(A®n>l>) l ) = (-) L+l '+ t ''^\ ' ? *' I V(iVA| | A| \N'0) (((N'O ® n"l") v > ® L) v \ \n'V 




where the harmonic oscillator matrix element is given by EQNS. (|T9")) - (f5!?)) and the final term is given by EQN. (139p . 
It is easily checked that if the initial state is a conventional meson then the expressions (|84|) - (|86)) collapse to the 
results of the previous section (|2"6l) - (|2"8]) . 

The amplitudes so obtained encompass many modes not yet discussed in the literature: these are too numerous 
to record here in their final polynomial form and so will be collated and discussed elsewhere. These hybrid decay 
amplitudes have the same generic functional form (|44]1 as the conventional meson amplitudes. The decays 

mP ± ^lP+ 1 S , (89) 

reproduce the results of REF. 0][23], where here and in what follows the label E for conventional mesons is omitted. 
The overall scale for these decays is set by the string overlap term; that scale factor is now being tested on the lattice 
and the agreement is good 0, Q. The authors of REF. Q found closed forms for decays of the negative parity states 
1IIP - but not for those with positive parity 11TP + . This arises from the mismatch between the orbital and radial 
characters of the wavefunctions, explicit in the wavefunctions (pB"|) and (jBT))) : the amplitudes derived here for HIP - 
states are finite polynomials in p whereas those for 1ITP + states arc a (rapidly decreasing) power series in p. 

The amplitudes derived here extend upon previous results in several different directions. The results are for 
arbitrary radial quantum number n' and so can be applied to more realistic radial wavefunctions that take account 
of the modified centrifugal barrier for hybrid states: any such wavefunction can be expressed as a linear combination 
of radial harmonic oscillator states. For the HIP 111 states modes such as 

mP ± -> lP+ 3 Si, and (90) 
mp ± -> lD+^o (91) 

have not previously been considered and have implications for hybrid width if hybrid mass is sufficient to allow decay 
Q; it may also be possible to test these on the lattice where masses can be adjusted to allow decay. The amplitudes 
derived here also include more exotic classes of decays that may be testable on the lattice, such as the decay of 
orbitally excited n hybrids or A hybrids: 

1IID ± — * E + E, (92) 
1AD ± -> E + E fee, (93) 

where here E stands generically for any conventional meson. Returning to the lightest HIP 111 states, the operator 
formalism admits a simple proof of the selection rule of Close and Page [|[ forbidding their decay to pairs of IS 
mesons, 

mP ± ^lS + lS. (94) 

Acting on a final state pair of identical IS mesons the decay operator is proportional to r, using EQN. (|5"3"|) : 

(((IS ® lS)o ® L)i||V||(A(gi ra'Z')j) = (((IS ® 1S) ® L) L \\V Tl + V r2 + V r/2 ||(A ® nT)i} (95) 

= (((lS«>lS)o®£)L||ri+r 2 +V r/2 ||(A<WZ');) (96) 
= (((lS®lS) ®L)i||r + V r/ 2||(A®n'0,) (97) 

Since neither the final state nor the operator have angular momentum about R, the initial state can only have zero 
angular momentum about R, corresponding to a E state: thus the decay of any A ^ state to a pair of IS mesons is 
forbidden: 

((IS ® lS)o ® L)i||V||(A ® n'l')i) =0 for A t^O. (98) 

The analogue of (|94[) can now be expressed in a more general form which is independent of the radial wavefunction 
and the total angular momentum I of the initial state: 

((IS ® 1S) ® £)i||V||(n® n'l')i) = 0, (99) 
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Since the zero appears for any n' it holds for any linear combination of n' such as a more realistic wavefunction which 
takes account of the modified centrifugal barrier for gluonic excitations. The result is also independent of the total 
orbital angular momentum I, forbidding the decay of not just I1P ± states but also their orbitally excited counterparts 
nD ± &c. As noted in REFS. @,[|, the rule is broken by outgoing mesons with different spatial wavefunctions: this 
can be understood by returning to physical units with different wavefunction widths /?i ^ P2, in which case the decay 
operator has a part which is proportion to ((3f — /3|)R and the decay from a IT state is allowed. A new result is that 
the rule also forbids such decays for hybrids with two phonons, 

((IS ® 1S) ® L)l||V||(A ® n'l')i) = 0, (100) 

and so forth. 

Isgur, Kokoski and Paton observed a more general rule forbidding the decay of UP^ hybrids to any pair of 
S-wave mesons with identical spatial wavefunctions, 

np ± ^s + s, (101) 

where here neither the initial nor final states are harmonic oscillator eigenstates. The present formalism admits a 
general proof of this rule, which turns out to be a special case of a more general rule: a state with A odd (even) is 
forbidden to decay to a pair of final states with the same spatial wavefunction coupled to In even (odd). The spatial 
wavefunction of any state |X) can be expressed as a linear combination of harmonic oscillator basis states 

\X)=J2^\n t k) (102) 

i 

The spatial part of the decay amplitude to any identical pair of states |X) is then 

(((X®X> 12 ®L)i f \\V\\{A®n'l')i) = ^c i c J -A y (103) 

where A^ is the amplitude (|84p with rijZj and njlj in place of n\li and 712/2- The terms Aij contain Talmi coefficients 

((^®n^), la |(JVA®n /, i ,/ )i' ia ) (104) 



which pick up (— )' + ' 12 under i j and are nonzero only for (— )'* +( i = (— ) A +' , using (jlip and (jT2|> . Together 
with the constraint (— = 1 it is easy to show that Aji = (— ) A+ll2 Aij and with the expansion p02[) above the 
selection rule follows: 

(((X®X)i la ®L) I/ ||V||(A®n / O ! ) = if (-) il2+A = -l. (105) 
So the decay of any II-state to a pair of S-wave mesons is forbidden, as they must be coupled to Z12 = 

(((X®X) il2=0 , 2) ...®L)/ / ||V||(n®n , / / ) i ) = 0. (106) 
A new result is that the decay of any conventional meson to an identical pair coupled to 1 12 — 1 is forbidden, 

(((X®X) Iia=1A „.®L), / ||V||(S®n / i / ),> = 0. (107) 

Few decay modes are entirely at the mercy of this rule since the recoupling (|23|) maps a final state of definite J12 onto 
a linear combination of Z12. The rule does, however, have implications for hybrid production, and so will be drawn 
upon in subsequent sections. 

VI. HYBRID MESON PRODUCTION 

Consider the most general case in which the final state can also include gluonic excitations with nonzero Ai , A2 : 

A^Ai+A 2 . (108) 

Recall that the flux tube quantum numbers Ai , A2 correspond to operators in the vectors r2 , i"i respectively. Acting 
with these operators on the final states {n! 2 l' 2 \, {n[l[\ gives a new pair of harmonic oscillators (n 2 'Z 2 '|, (n'{l"\. In this way, 
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any amplitude involving gluonic excitations in the final state can be expressed as a linear combination of amplitudes 
involving only conventional mesons: 

((((Ai ® n\l' 1 ) h ® (A a ® n^ahJiu ® L )*/II V II( A ® »'0i) = ® "i^'i ® ( Aa ® "a^iJlulK'ii ® <^)h 2 ) 

(((«/'/ ® n^')( 12 ) ® £)i,||V||(A ® n'/'Jj), (109) 

where the second term is that calculated in the previous section and summarised in (|84|) - (|86p . The requisite trans- 
formation is found by recoupling the Ai, A2 operators so as to interchange their positions, and using 

KASn'iV,*) = \n"l")^-{n"l"\\K\\n'l'). (110) 
lb" 

For the general case the recoupling is 

(((Ai ® nil£)ii ® (A2®n' 2 l' 2 ) h )iJ(n'{l'{®nZl%) ll2 ) 

' l[ Ai h 

( _ ) A 1+ A 2+ i 1 n hi2 (_y" j a 2 i> i 2 Wn'/r/iiA^in^DK^iiAiiK;^). (in) 



l 2 

ill pi 1 

If one of the final state mesons is a conventional meson, 
(((Ax ® nili)j, ® (S ® n' 2 !2)( 2 )iiJK''i ® ^a'k.) 

= (-)' 1+/l2+ ' 2 n ;i j £ j» j (n^a'HAilK/O^^ir/, (112) 

and if both final state mesons are conventional mesons the recoupling is trivial: 

(((£ ® n'J'^ ® (E ® n^) la )i J("i7i ® «£#)««) = fyn'/fyri'^W^i'i'- (113) 

It follows immediately that these more general amplitudes recover the results of the previous section if both outgoing 
mesons are of the conventional type. 

The results so obtained are the first calculation of hybrid production amplitudes. Processes such as 

E -> n + E, and (114) 

e -► n + n (115) 

will be discussed elsewhere in the context of e + e~ i?-factories Q. The next section considers a special case of 
such processes in which a 1IIP - hybrid is produced along with a IS state. In this case the recoupling (|112|) is 
straightforward, mapping onto a pair of IP states coupled to £12 = 1: 

((((n ® lP)x ® (E ® 1S) )i ® L)i,||V||(A ® n'1%) oc (((IP ® 1P) X ® L) lf ||V||(A ® n'J'Jj), (116) 

so that, using (|105[) . a new selection rule arises 

E^mp- + is. (117) 

Another way of seeing this is to replace IT by the operator r 2 and for the harmonic oscillator in ri that | IP) oc ri|lS), 
thus 

((((n ® 1P)i ® (E ® lS) )i ® L)i f \\V\\(A ® n'l')i) 

cx ((((r 2 ®n)i® (IS® 1S) )i ® i)i / ||V ri +V, 3 + V r/a ||(A(8)n , 0i) (118) 
The decay operator then gives two terms, one proportional to 

(((R ® (IS ® 1S) )i ® L) lf 1 1 (A® n'l')i) (119) 

and another proportional to 

((((R ® r)i ® (IS ® 1S) )i ® L) h ||r + V r/2 ||(A ® n'l')i), (120) 
both of which arc zero for an initial E state on the account of the vector R in the final state. 
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FIG. 4: Some possible sources of hybrid states in e + e B-factories. 

VII. THE SUPPRESSED PRODUCTION OF THE 1+ HYBRID 

There has been speculation that the surprisingly large double cc cross sections in e + e~ £?-factories could prove 
a copious source of charmonia hybrids, or likewise the decay of physical bb states. Analogously, it is possible that 
light quark hybrids could be buried in the wealth of data anticipated on charmonia decays. Such processes will be 
considered in more detail in REF. Q; in this section it will be shown that within the framework developed there the 
production of the 1HP _ hybrid along with a conventional IS meson is forbidden. 

The processes under consideration share common features as depicted in FIG. [U annihilation into a virtual photon 
or gluons and the creation of a virtual quark pair which then decays via pair creation. In the approach of REF. 
only the second part of the process is calculated - the strong decay of a virtual quark-antiquark pair into the final 
state pair. In a non-relativistic potential picture, the pair that emerge from the black box can be any admixture of 
orbital, radial or gluonic excitations consistent with the initial state quantum numbers, and a priori the heirachy 
of mixing angles between such states is unknown. Arguments given in REF. Q suggest that gluonically excited 
modes are suppressed relative to E modes. As an example consider the 1 channel, whose intermediate states could 
be admixtures of |n 3 ESi), |n 3 EDi) and |n 1 IIP^). Any significant II admixture would be in conflict with known 
prevalence of experimental modes such as 

e+e" -> rj c ip, (121) 

V -> P7T, (122) 

V> 6itt, (123) 

the first two on account of the selection rule Q101[) and last because of the spin-singlet selection rule, a zero in the 
9-j coefficient in (|2~5j) for s = s\ = s 2 = 0. This suppression may be expected because inside the black box of FIG. Q] 
will be factors sensitive to the wavefunction at the origin, suggesting |n 3 ESi) states will dominate. For channels with 
quantum numbers other than 1 and ^ both the E and glue-excited states have nodes at the origin so equivalent 
arguments do not apply, but there is phenomcnological evidence suggesting the E states dominate. For \i decay, for 
instance, the intermediate state could be populated by |n 3 EP2) and |n 1 IIDj) but the prevalence of modes such as 

\2 -» <H>, UUJ, 7T7T, 7777 (124) 

would argue against the importance of the latter, for here the selection rule (|106[) is expected to be exact. Deferring 
further arguments to REF. [3], suppose that the gluonically-cxcited modes can be neglected, in which case the spatial 
wavefunction of the pair emerging from the black box can be any admixture of E states consistent with the quantum 
numbers of the initial state. All such processes as depicted in FIG. [4] are then be subject to the selection rule (|117|l . 

E^mp- + is, (125) 

valid in the limit of equal width wavcfunctions for the final states but independent of the radial wavefunction of 
the intial state. This has implications for the lightest hybrid states |1 1 IIP^) with 1 quantum numbers and the 
|l 3 nPo" 12 ) with (0,1,2) h quantum numbers. In charmonia decay, for instance, some channels whose quantum 
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numbers arc ripe for the production of the exotic 1 + should be forbidden in this limit: 

tp -f* nip in P-, F-wave, and (126) 
Xi 7Ti7r in S-, D-wave. (127) 

Strong observation of such modes in the ~k\ (1400), ~k\ (1600) channels would argue against a hybrid interpretation for 
them, although the different spatial wavefunctions of the 7r and ~k\ may lead to significant breaking of the selection 
rule |7|. For charmonia states the rule would be expected to hold more strongly: this would argue against a hybrid 
interpretation for the X(3940) observed in 

e+e- -> J/^X(3940), (128) 

if it is found to have J ^ quantum numbers, but not neccessarily if it has 1 ++ . A hybrid interpretation is thought to 
be untenable in any case, on account of the low mass. On the other hand if the 1" y(4260) is a 1 1 1 nP 1 ^ ) hybrid then 
it should have (0, 1, 2) h partners |l 3 IIPg 1 2 ) in the same mass region; because of the selection rule, the absence of 
such peaks in the data is understandable. The selection rule may help discriminate between hybrid and other exotic 
or non-exotic interpretations of heavy or light-quark states, include those with non-exotic quantum numbers. For 
instance, searching for forbidden modes such as 

T,e+e- y(4260)ry c (129) 
may help discriminate intepretations of the F(4260). 



VIII. CONCLUSION 



The Talmi-Moshinsky transformation has been applied to the strong decay problem and several new results derived. 
The method allows easy derivation of previous results for conventional mesons and extends upon those to encompass 
decays involving any orbital and radial excitations. By factorising the string degrees of freedom amplitudes for decays 
involving gluonic excitations are a simple generalisation of those involving only conventional mesons. Many of the 
amplitudes so obtained have not previously been considered; their final polynomial forms will be collated elsewhere 
and applied to conventional and hybrid meson production in charmonia decays Q ■ 

Within this formalism the well-known selection rule against the decay of a hybrid to a pair of identical IS mesons 
has a natural explanation. Using the Talmi-Moshinsky transformation it is possible to generalise this rule to the 
case of identical wavefunction that need not be harmonic oscillators, in agreement with numerical calculations with 
realistic wavefunctions. There appears a new class of selection rules which have implications for the production of 
1 h exotic hybrids and may help discriminate interpretations of states with non-exotic quantum numbers. Owing 
to its being purely geometrical in nature, the Talmi-Moshinksy transformation may find application to other quark 
model processes with the geometry of FIG. There are generalisations of the Talmi-Moshinksy transformation for 
systems described by more than two Jacobi coordinates and these may be useful in, for instance, meson rescattering 
amplitudes or processes involving baryons. 

The flux-tube model naturally suggests a decay operator with pair creation in spin 1 and provides a parameter-free 
scale for processes involving hybrids relative to those involving only conventional mesons. Lattice QCD is now able 
to calculate strong decays and these two features of the model appear to be confirmed [|| . Comparing to the lattice 
allows models to be tested and refined; at the same time models can shed light on the dynamics of the underlying 
QCD of the lattice. Further comparisons of the type [f| can be made easily using the results of the present work. 
The formalism developed here could also be applied to decay amplitudes via the 3 Si operator and to more general 
processes allowing for spin- or orbital-flips at the pair creation vertex. 

Generalisations of the operator approach employed here may be of use to the problem of rescattering between hybrid 
and conventional meson states, using the ideas of [l9| . Likewise it may be possible to employ similar techniques to 
model flux-tube de-excitation decay modes, predicted in lattice calculations to be significant for hybrid states [2fj| : 
this could help in the interpretation of the Y(4260) state. 

The recoupling Q109p that relates the spatial wavefunction of outgoing hybrid states to those of conventional 
mesons can also be performed in the \(s®nl)j) basis. In this way, final states involving hybrids can be expressed 
as a superposition of conventional meson states and their production amplitudes relative to conventional mesons is 
independent of the initial state; this is applied to hybrid pair production elsewhere Q ■ 

That the wavefunctions in one set of Jacobi coordinates translate into a finite number of wavefunctions in the other 
is a feature peculiar to harmonic oscillator wavefunctions. The calculations presented here could be used as a basis for 
calculations with more realistic wavefunctions: the coefficients in the translation from these to the harmonic oscillator 
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basis are well-known. It is straightforward to generalise these calculations to allow for different wavefunction width 
in the initial state or to include a localised pair creation region with spherical symmetry: these would appear only in 
the final integrals (|4T|) and (|42|) . Those integrals could also be performed numerically to incorporate the cigar-shaped 
pair creation region of [l3j . 
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